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1.  Introduction 


The  classical  theory  of  ignition  in  energetic  substances  is  based  on  the  appearance  of  hot  spots. 
Hot  spots  are  the  localized  areas  with  temperatures  much  higher  than  the  average.  These  spots 
can  be  triggered  by  different  physical  effects  and  processes.  One  of  them  is  based  on  the 
mechanism  of  adiabatic  loading.  The  classical  theory  of  this  mechanism  for  the  liquid/gas  media 
is  presented  in  the  monographs  in  references  1-4.  The  main  element  of  this  mechanism  is  that 
the  easily  compressible  gaseous  bubbles  under  given  pressure  become  hotter  than  the  less 
compressible  condensed  liquids  surrounding  them.  In  the  first  paper  of  this  study  (5),  we 
demonstrated  that  the  classical  theory  should  be  essentially  revised  by  considering  the  effects  of 
the  surface  tension.  The  revision  mostly  focuses  on  the  influence  of  the  dissolved  bubbles  of 
small  sizes. 

In  the  second  paper  of  this  series  (6),  we  explored  spherical  holes  in  solid  materials.  These  holes 
can  serve  as  the  local  hot  spots.  They  are  completely  different  mechanisms  than  the  local 
temperature  concentrator.  The  gas  inside  the  hole  plays  no  role  in  this  mechanism.  The  key  idea 
associated  with  this  mechanism  is  as  follows.  When  the  flammable  liquid  is  subjected  to 
external  impact,  the  pressure  is  totally  transferred  from  the  liquid  to  the  dissolved  bubble.  The 
situation  changes  completely  when  the  liquid  is  substituted  with  solid.  The  solid  then  plays  the 
“shielding”  role  and,  at  least  partially,  precludes  direct  transfer  of  the  external  impact.  Thus,  the 
local  adiabatic  stressing  of  the  solid  remains  the  only  source  of  the  local  temperature  rise. 

One  key  difference  between  liquid  and  solid  energetic  materials  was  emphasized  in  Grinfeld  and 
Bjerke  (6).  Namely,  at  equilibrium,  the  gaseous  bubble  in  the  liquid  always  has  the  spherical 
shape.  At  the  same  time,  it  can  and  does  have  arbitrary  shapes  in  solid  substances.  Thus,  the 
analysis  of  ref  6  is  insufficient.  Therefore,  in  this  third  part,  we  explore  the  role  of  the  form 
factor  of  the  hole.  In  order  to  make  a  more  feasible  and  transparent  study,  we  limit  ourselves 
with  the  two-dimensional  case.  We  also  dwell  on  the  explicit  analysis  of  the  elliptic  holes. 

Our  analysis  reveals  one  element  of  key  importance — the  “shielding  effecf’  disappears  for  the 
ellipses  of  small  eccentricity.  Thus,  intensive  heating  of  gases  becomes  important  again  to  fill 
these  holes.  We  establish  explicit  relevant  equations/formulas.  Also,  we  present  the  exact 
nonlinear  formulation  of  the  problem.  The  nonlinear  system  can  be  solved  numerically.  For 
those  numerical  solutions,  the  established  exact  solutions  can  be  used  for  verification  purposes. 
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2.  Statement  of  the  Problem 


Exact  Nonlinear  Statement  of  the  Problem 

In  figure  1 ,  an  isotropic  matrix  contains  an  elliptic  inclusion,  which  usually  plays  the  role  of 
stress  concentrator. 


Figure  1 .  An  elastic  plane  with  a  hole. 


Our  notation  follows  the  monographs  in  references  7  and  8.  We  shall  use  the  notation  x(x)  and 
X  (x)  for  the  radius  vectors  of  a  material  particle  with  Lagrangian  coordinates  x'  in  the  initial 
configuration  and  actual  configurations,  respectively.  The  radius  vector  X (x)  can  be  presented 
in  the  following  form; 

X  =  x  +  w(x) ,  (1) 

where  m(x)  is  the  displacement  vector. 

Then,  we  use  x.  (x)  =  Sx(x)/5x'  and  X.  (x)  =  x.  +  dul dx‘  for  the  covariant  bases  in  the  initial 

and  reference  configurations.  The  notations  x.j  =  x.  •  Xj  and  X^  =  X.  ■  X .  are  used  for  the 

covariant  metric  tensors  of  the  initial  and  reference  configurations,  respectively.  For  the 
contravariant  metrics,  we  use  x^  and  X'^ ,  respectively.  We  use  T'  and  T.  p  for  the  covariant 

differentiations  of  tensors  in  the  initial  and  reference  configurations.  Tensors  x.j ,  x'^  are  used  for 

lowering  and  raising  indexes. 

We  use  u'  (x),  w,  (x)  for  the  components  of  displacements  with  respect  to  the  basis  in  the  initial 
configuration,  i.e.. 
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(2) 


u  =  ux^  (x)  =  u^x'  (x) . 


Note  the  following  geometrie  relationships: 

Xj  =  a  'jX.  =  {S‘  +  Uy)x, ,  Xj  =  b'jX.,  =  d‘,  , 

and 

Xij  =  X,J  +  2Uy  ,  U,J  =  +  M^.|,  +  ) , 


(3) 

(4) 


where  m,,.  is  the  tensor  of  finite  differenees. 

We  use  Po  and  p  for  the  mass  densities  in  the  initial  and  aetual  eonfigurations.  In  addition  to 

their  physieal  meaning,  they  make  eonneetions  between  the  divergenees  of  any  eontravariant 
tensor  T‘  in  the  initial  and  aetual  eonfigurations: 

r 

Po 


Pq  'j-’i 


(5) 


We  use  P‘^  and  p^‘  for  the  Cauehy-Green  and  Piola-Kirehhoff  stress  tensors,  respeetively.  These 
two  tensors  are  eonneeted  by  the  following  relationship: 

^  p^ipi,  (<5,  p,  ^  P_J,„  (  )x-‘ .  (6) 

P  Po 

The  exaet  universal  bulk  equilibrium  equation  within  the  matrix  is 

=  0  .  (7) 

The  exaet  universal  equilibrium  equation  at  the  interfaee  matrix/gas  is 

—  P'"  .  (8) 

Pq 

When  the  matrix  deforms  elastieally,  the  substanee  is  fully  eharaeterized  by  the  internal  energy 
density  e{u^^j,ri^  per  unit  volume  in  the  initial  eonfiguration;  rj  is  the  entropy  density  per  unit 

volume  in  the  initial  eonfiguration. 

We  assume  that  initial  distribution  of  speeifie  entropy  rj  is  uniform — //(x)  =  77^^,  =  const — and 

remains  unehanged  under  the  loading.  We  make  the  same  assumption  about  the  gas  or  liquid 
within  the  hole:  //(x)  =  77^^,  =  const .  In  this  ease,  exaet  equilibrium  equations  7  and  8  ean  be 

rewritten  in  the  following  form: 
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(9) 


5m, I, 


=  0. 


J\i 


and 


P_ 

Po 


du 


q\] 


-Pgas^k  • 


We  assume  that  displaeements  U‘  are  speeified  at  the  external  boundary  of  the 
heterogeneous  solid; 

m‘|  =U'. 

'S, 


(10) 


(11) 


To  make  equations  9-1 1  mathematieally  elosed,  we  also  need  to  know  the  dependenee  of 

P gas  ~  P gas^ gas  gas^  • 


3.  Semi-linear  System  of  Adiabatic  Loading 


The  system  just  presented  is  exaet;  however,  it  does  not  allow  deeper  mathematieal  analysis.  It 
will  be  analyzed  further  by  eomputer-based  modeling.  Computer  modeling,  in  turn,  requires 
eomparison  with  some  exaet  solutions  and  a  qualitative  understanding  of  the  studied 
phenomenon.  Therefore,  we  proeeed  with  some  simplifying  assumptions.  Namely,  we 
eoneentrate  on  the  phenomena  for  whieh  the  approximation  of  linear  elastieity  provides  the 
physieally  adequate  deseription  of  the  matrix’s  behavior. 

In  the  approximation  of  linear  elastieity,  equations  9  and  10  ean  be  replaeed  with  the  following: 

(12) 


and 

C^V,M.=-n(F  ,11  ), 

r}  kj  j  ^  \  gas^  I  gas  J  ’ 


where  is  the  tensor  of  adiabatie  modules. 


(13) 


Adiabatie  loading  is  aeeompanied  by  the  ehange  T^^=6-  6^  of  the  absolute  temperature  6^ , 
whieh  is  deseribed  by  the  following  approximate  relationship; 


ad 


'Y  _  ^0  Qkl 

~c 


k,l  ’ 


(14) 
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where  C„  is  the  speeifie  heat  eapaeity  at  fixed  deformations  and  B'‘‘  is  the  tensor  of  eoeffioients 
of  thermal  expansion. 

The  actual  volume  of  the  gas-lilled  hole  will  be  approximated  by  the  following  relation; 

(15) 

s 


The  case  of  isotropic  matrix  is  characterized  by  the  special  structure  of  the  tensors  and 
B^‘: 


A 


V  '7 
B‘^  =  -K,ax^ , 


2 

l-2v 


2(^7 


(16) 


where  and  are  the  Lame  adiabatic  elastic  modules  ( Xg  and  iUg=  iu^=  ju  are  the  Lame 
isothermal  elastic  modules),  Kg=Xg  +  2jUgl3  is  the  isothermal  module  of  volumetric 
compression,  and  a  is  the  coefficient  of  thermal  expansion. 

For  isotropic  substances,  we  arrive  at  the  following  relationships: 

=  (17) 


and 


■ciP''  =XgX'^u';,+2uulf,, 


the  second  of  which  can  also  be  presented  in  the  following  form; 


adP"  = 


Xg  + 


c 


+2jUgM^y 


(18) 


(19) 


4.  An  Elliptic  Inclusion  in  a  Stressed  Isotropic  Plane 


In  this  section,  we  focus  on  one  rich  class  of  problems — permitting  explicit  solutions.  These 
problems  deal  with  an  elliptic  hole  within  an  infinite  isotropic  matrix,  with  specified 
displacement  gradients  k.j  at  infinity. 

u.  KyX^  at  I  x“  I  ^  00  .  (20) 
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We  will  be  looking  for  solutions  of  the  system  (equation  20)  next. 
Within  the  matrix: 


U^  =  K..X^ - (j)  ■  + - ^ijk 


Within  the  inelusion: 


u^  =  a)^■T 


where  is  a  certain  symmetric  matrix  and  ^zi(x)  and  <^(x)  are  the  harmonic  and  bi-harmonic 
potentials  of  an  ellipse  having  the  shape  of  the  hole: 


(*  +  1  r  *  1  ^  * 

^zi(x)=  da  In  — ,^(x)=  da  —  In  — +  1  ,r(x,x)=| 

J  0)  y  J  CO  2  '  ^ 


x-x  . 


The  potential  ^(x)  and  its  first  derivatives  are  continuous  across  the  ellipse  boundary.  The  same 
is  true  about  the  potential  ^(x)  and  its  derivatives  up  to  third  order.  The  following  relationships 
for  discontinuities  of  the  higher  derivatives  are  valid: 


Inside  the  ellipse,  the  second  derivatives  of  the  potential  (j)  are  constant  and  given  by  the 
following  formula: 

where  is  the  length  of  the  semi-axes  of  the  ellipse  and  is  their  director. 

The  analogous  formula  for  the  fourth  derivatives  of  the  potential  ^(x)  inside  the  ellipse  is  as 
follows: 


iTva^ia+laA  2na,a.  InaAa.+laA 

K  =-M  - _ Tkd _ 11]  I  I  I _ 12  T _ 1 V  2 _ \_!_J  J  J  J  ^261 

h,ijkl  ^"ijkl  ,  .2  'irii'u'l/  ,  \1  ijkl  /  x2  '2!'27'2/t'2Z  ’  ) 

(aj  +  aj  j  (^1  +  <^2  j  +  ^2  j 

where  the  tensor  is  equal  to 

^ijki  ~  hih Alkali  hJijhk^ii  hAAikhi  +  hihjhkhi  +  ^nhAikhi  +  ^lAjhkhi  •  (22) 

Direct  verification  gives  us 
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(28) 


= 


Ina^  (flj  +2a2) 


.  x2  '■Ik'-ll  ,  .. 

(aj+aj) 

/  +±^l  I  =1N 

nkhl^  .  kr 


4  Z.4  7 


Ikum. 


1^*2 


\hkhl 


2;raj  (aj  +2aj) 
(flj  +02) 


2  ^2i:^2Z 


aj  +a2 


aj  +a2 


Combining  equations  21  and  24-26,  we  get  the  following  formula  of  the  displaeement’s 
gradients  on  the  ellipse’s  boundary; 


'•ds 


1 


1  -  1 


A 


^.i  ‘Pjk  +; 


E 

h.iiki 


=7^.+  — r;-Ad,- 


8;r//(l-K)  ■■ 
1 


/s+o 

1 


1 


Inn  '  8;r//(l-v')  ^  2//(1-k) 

Using  equation  29,  we  arrive  at  the  following  formula  of  the  stresses  in  the  matrix; 


(29) 


mn  mn  , 

Ps  =Po,  + 


i - 1 - ^nrn  ^  ^in  .rn) 

2(l-u)  n  nA{\-v\  "  n 


^kl 

T  T  fl  flh. 


1  (  m  n  , .  c*mn  \ 

Y~^\^  n  -v^S  )n,n,T 


kl  _kn  m 

—  X  n  n 


k  ’ 


(30) 


where 


mn 

P^  =P 


2k 


vl-2^ 


rj  ij  mn  ,  im  nj  ,  in  jm 

- — X  X  +X  X  -\-X  X 


K,.. 


(31) 


Contraeting  equation  30  with  ,  we  get 


^mn^^  ( _ mn  .  Anm..^kl\^^ 

Ps  ««  =i^co  +^..kl 


(32) 


where 


Ainn..  ^ ri  mn  ^  at"  ^  ^  urntn  ,  ^  Alin  c*m)  c*n  c*m 

Ak,  =  — - -Nk, - — - -M  +-N\S,  ’-S,S,  . 

2(1 -u)  n  n4(\-vj  n 


(33) 


In  view  of  the  boundary  eondition  (equation  13)  and  the  arbitrariness  of  the  normal  orientation 
that  equation  32  implies, 

^mny  ^  ^mn  ^  _  ^4) 

Equation  34  should  be  treated  as  a  linear  equation  with  respeet  to  the  unknown  tensor  A . 

Taking  the  skew-symmetrie  part  of  equation  34  and  using  equation  33,  we  get 
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(35) 


\ki]  =  0  • 


Thus,  the  tensor  is  necessarily  symmetric.  For  the  symmetric  part  of  tensor 
equation  34  obviously  implies  the  following: 

=-pT  -  p{V)S"'\ 

For  the  actual  volume  V ,  we  get 

F  =  F,  l  +  < - ^ . 

[  2;tp  ’  8;r//(l-v^) 

Relationship  37  is  implied  by  the  following  chain: 


F  =  Fn  +  f  dSu^.n  =  Fn  +  f  - - - 

^  °  J5  I  "  27:p  %7:p{\- 


■■V„+ldS 


%7rp{\-v„) 


im  r  n 


=  Fo  +  f  dF\K\ - ^r^'‘>,y+ - - 


=  f;  l  +  < - ^ . 

iTvp  ■■  8;z-//(l-v^,)  ■■  " 

Using  formulas  24,  and  28,  we  can  rewrite  equation  37  as  follows: 

F  1-2U 

-=i+/^:+T7^ — , 


A[\-v^)7rp 


as  it  follows  from  the  following  chain: 


F  .  1  ..  1  ... 

—  =  \  +  < - ^‘dk 

K  '  2;z-//  ■■  ’  87rp(l-v^)  ' 

=  1  +  /f i  +  —  N,, - i =  1  +  <■  +  N, . 

2np  ^np{\-y^)  ^{\-v^np 


The  following  explicit  formula  of  derived  in  appendix  A: 
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1 


Amn)..  _ 

A(kl)  - 


4a^a 


2 _ /(»!?«)?  ? 

^  tl  t')  Ia/kI' 


2  n  ^2  n(/t^2/) 


2(1“^^)  («!  +a2) 

_  2(1  -  )«!«!+ (l  -  2k J  ajaj 

2(l-^;,)(«i+«2f 

-  (l  -  212^  ]  a^a. 


I  I  /^/”  + 


Iva^a^  “  (l  “  '^^rj ) 


2(1-^^)(«i+«2) 


/  /  /'”/” 
2  nknri  ^2 


2(1-^^)(«i+«2) 


ujjTi’:  - 


2  ^2/t^2/‘'l  n 


ajaj  (l  -  2k^  ]  +  ) 

2(l-'^^)(«i+«2f 


/2./2//2"/2-  (41) 


Let  B  be  a  tensor  symmetric  in  the  upper  and  lower  indexes  such  that 

x(mn)..  nkl..  —  S2(m  cn) 

^..(kl)  ^..pq  -^P  ■ 

Tensor  is  calculated  explicitly  in  appendix  A  and  is  given  by  equations  A-3  and  A-4. 


(42) 


Then,  equation  36  implies 

t'"'=-B'-'JipZ-+p(V)S-].  (43) 

Using  equation  33,  we  can  rewrite  equation  30  as  follows: 

mn  mn  ,  (  Amn..  ,  <pn  <pm\  kl  km  n  ,  1  /  m  h  <pmn\  kl  kn  m  i  A  A\ 

^4  +[Aki  +^k^i  nnk+— — [n  n  -v^d  jn^riiT  -r  n  .  (44) 

At  last,  using  equation  36,  we  can  transform  equation  44  as 


n  n  -v^o  ]n,^n,T  —r  n  n,^ 


„mn  „mn  ,  (  Amn..  ,  i?«  i?ffi  \  ^  £mn\ 

Ps  =Poo  +[A..ki  +^k^i  jU  nf^k+Y^[f^  n  -k/  j 

^'7 

/Ar\  ^mn  ,  nm  aa.  k(m  n)  ,  1  (  m  n  ^mn  \  kl 

=  -p{V)d  +r  -2r  n  rij^A - [n  n  -v  o  }ni^n,T  . 

l_y^\ 

Let  us  consider  any  vector  tangential  to  the  unit  normal  n. .  Then,  equation  45  gives  us 


(45) 


p:xtn=-piv)+T: 


1 1  — 

m  n 


1-K„ 


(46) 


Using  equations  A-3  and  A-4  from  appendix  A,  we  can  rewrite  equation  49  as  follows: 


77||7i*2:/,p/|^  +5„/f/i72^4^  +7321/2/:/, + 


[a"  +  p(V)S«'_. 


(47) 


or  else  as 
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(48) 


+p)-  Bjx  + p) 

-  Bj‘j,  (<7r/„;„ +p)-  Bj‘A  (<Tr4,4, + p)-bi;‘i'Aj. 


2?)^«  • 


Equation  39  now  takes  the  following  form: 


v-K_  1 


3  ,  3 

-\-a2 


Vg  ju  a^a2{a^+ a^)' 


v^oy 


.  l-2u 


+  P.i. 


2// 


2{\-v^^a\  +{\-2v^^{^a^al  - af aj )  + 

2jj.a^aj^  (aj  +a2) 

2  (l  -  +  (l  -  2v^  j  -  a\a^  j  +  2Ka2 


pi%u 


2iia^aj^  i^a^+a^) 

Introdueing  the  eeeentrieity  a^l =  ^  ,  we  ean  rewrite  equation  49  as  follows: 


V-V,  1  1-s  +  s 


2  ry\ 


K 


P 


v^oy 


=  P.k 


,  ('-%)«’ +  ^;>r 


2// 


-  +  - 


•(!  +  .) 


P 


hph. 


l-K„+l/„f%^^  .  .  (\-2v^){\-s)  pPJ 


(49) 


s{\  +  s^  /u 


■4^4? +■ 


l  +  e 


2p 


{hpkq  hphq)-  (50) 


5.  Some  Special  Cases 


In  the  ease  of  a  eireular  inelusion,  eoeentrieity  s  is  equal  to  1  and  equation  50  reads 


v-v,  1 

V,  p 


^v)  1-E 


V^oy 


P 


r  oo 


{juhj  •- 


(51) 


or  else 


v-v._^ 

K  p 


V^oy 


+ - (^l«+^2«) 

M 


(52) 


where  cTj^  and  cr2„  are  the  prineipal  stresses  at  infinity. 


^lo.=Plk.kj^^lo.=Plhkj- 


(53) 


At  =  ^  «  1 ,  we  ean  reduee  equation  B-3  in  appendix  B  to  the  following  form: 
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V-Vg_l  f  v]_(i-2v,)£  +  v,  (1-2k,)^  +  1-i/,  cr. 


.  (54) 

2// 


In  the  lowest  order  terms  in  s  ,we  get 


P  —  =-y^ioo-(l-^)^2co- 


5.1  The  “Linear”  Gas  Model 

For  further  insight,  let  us  eonsider  the  simplest  case  of  a  “linear  gas”  defined  by  the  following 
relationship: 

V  -V 

P  =  K^^.  (56) 

When  dealing  with  the  circular  hole  and  the  model  of  linear  gas,  equation  52  gives  the  following: 


v,-v 

K  m  +  k,  ■ 


(58) 

M  +  K, 

where  is  the  pressure  at  infinity,  defined  as  =  -(cj^  +  . 

Usually,  the  adiabatic  incompressibility  of  gas  is  much  smaller  than  the  shear  modulus  p  of 

the  matrix.  Then,  equation  58  teaches  us  that  there  is  a  strong  “shielding”  factor  O  : 

<b  =  2Kp-v^)l(^  +  K,).  (59) 

In  other  words,  the  induced  pressure  p  of  the  gas  is  much  smaller  than  the  pressure  p^  acting  at 
infinity,  hence,  the  smaller  induced  temperature. 

For  the  case  of  the  arbitrary  elliptic  hole  filled  with  the  linear  gas,  equation  50  reads 


l-2v/„ 


l-2v/„ 


where 
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(61) 


K  (\-S+£"]l+S 

r  =  ^,  A  =  -i - 

JU  \  +  £ 


’-y 


For  finite  F  and  ^  ^  0 ,  equations  60  and  61  give  us  the  following  in  the  lowest  order  terms: 

;?  =  -^,^ioo-(l-^)^2oo-  (62) 

Equation  62  shows  that  for  the  eraek-like  hole,  the  shielding  effeet  of  the  matrix  disappears 
completely;  the  pressure  with  the  hole  is  the  same  order  of  magnitude  as  the  stresses  applied  at 
infinity. 

5.2  The  Model  of  Ideal  Gas 

Equation  55  shows  that  the  formula  (equation  62)  for  the  gas  pressure  is  valid  not  only  for  the 
model  of  the  “linear”  gas  within  crack-like  holes  but,  actually,  for  any  gaseous  liquid  media. 

We  are  reminded  that  the  adiabatic  processes  of  the  ideal  gases  are  described  by  the  following 
relationships: 


(63) 


y 

e 

A 

W) 

9 

iv) 

The  constant  /  is  defined  as  y  =  /  C^, ,  where  and  Q  are  the  specific  heat  capacities  at 

fixed  pressure  and  volume,  respectively. 

Thus,  for  this  model,  equation  50  can  be  rewritten  as  follows: 


=\-p  ,  (1-^)^’  +  '-',  <T, 


{- 

Y-\ 

-E±\. 

-£^£^ 

f-1 

J 

£ 

1  -  +  V£^ 


£{\  +  £^  /U 


£{\  +  £)  /U 

where  is  the  original  pressure  of  gas  in  the  hole. 


+ 


\  +  £ 


Ip 


(64) 


6.  Discussion  and  Conclusion 


We  formulated  exact  nonlinear  systems  (equations  9-1 1)  to  numerically  explore  the  stress- 
induced  hot  spots  around  holes  of  different  shapes  filled  with  gases  or  liquids. 
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We  then  formulated  semi-linear  systems  (equations  12-15)  based  on  the  assumption  of  small 
deformations  in  the  solid  matrix.  The  “pressure-volume”  relationship  for  the  gas  or  liquid 
remains  nonlinear. 

With  the  help  of  the  semi-linear  system,  we  analytically  and  explicitly  explored  the  problem  of 
thermoelastic  stresses  and  deformation  in  an  unbounded  isotropic  matrix  weakened  by  the  gas- 
or  liquid-filled  inhomogeneity. 

The  analysis  showed  two  opposite  effects.  First,  in  the  case  of  a  circular  inclusion,  there  was  the 
“shielding”  effect.  In  other  words,  due  to  the  large  shear  resistance  of  the  matrix,  the  pressure  of 
the  gas  was  considerably  smaller  than  the  stresses  in  the  matrix.  This  effect  does  not  exist  when 
dealing  with  liquid  energetic  materials.  Therefore,  contrary  to  the  case  of  liquid  energetic 
substances,  when  dealing  with  solids,  the  spherical  holes  had  a  much  less  chance  of  becoming 
hot  spots  (at  least,  for  relatively  small  stresses). 

In  the  crack-like  holes,  the  “shielding”  effects  from  shear  resistance  disappeared.  The  gas 
pressure  within  the  hole  was  the  same  order  of  magnitude  as  the  stresses  in  the  matrix. 

Therefore,  the  behavior  of  gases  within  these  holes  was  completely  analogous  to  the  behavior  of 
gases  within  spherical  bubbles  in  liquids.  The  gas-filled  holes  in  solid  explosive  materials  had 
better  chances  of  becoming  the  “hot”  spots. 
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Appendix  A.  A  Closer  Look  at  the  Tensor  A*-"'"},'  and  Its  Symmetric  Inverse 


Tensor  can  be  presented  in  the  following  form: 

j{mn)..  _  j  imini  J  i  A  ]mini  I  i  A  ]^]ni  J  i  A  ]mini  i  ]  ( A-]^ 

^  ..(kl)  ^iri  h  nkhl  ^  n  ^2k^2l  ^  ^2r2  ^2  nkhl  ^  ^22^2  ^2  ^2k^2l  ^  ^2  ^Uk^2l)  *  -*-/ 


The  set 


f-  — {imini  1  imini  j  /"«/«/  /  jmjnj  j  l{min)i  j  ) 

vj-  _  hkhi  5  h  hhkhi  ’  h  hhkhi  ->  H  h  H(khi){ 


(A-2) 


is  called  canonical  basis. 

Let  be  the  inverse  tensor  with  respect  to  A^lf"  in  the  sense  of  equation  42  from  the  body  of 
this  report.  We  assume  that  the  tensor  is  symmetric  with  respect  to  the  upper  and  lower 
indexes.  We  will  be  looking  for  in  the  same  form  as  follows: 

nkl..  _  D  J  -L  R  J^J^J  J  -I-  R  J^J^J  J  -I-  R  J^J^J  J  1  (A 

^  ..pq  ^\\^\  hhphq  ^12^\  h^2p^2q  ^2\h  hhp^lq  ^  22^  2  ^  2^  2  p^  2q  ^2  ^\{p^2q)  * 


Direct  calculations  lead  to  the  following  formula  of  the  coefficients  B^^  ,B^^  ,B^^  ,^22 ,  and  b  : 


5,,  =- 


l-K  2(l-K)a2 +(1-2k)< 


l-2v 


■’^2  ~ 


l-K  2Ka2  -(l-2i/)aj 


l-2v 


5.,  =- 


l-v  2va^  -(1-2k)< 


l-2v 


"’^22 


l-v  a2(l-2K)  +  2aj(l-K) 


l-2v 


h  =  2(1-1/) 


(flj  +  ^2)^ 


(A-4) 
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Intentionally  left  blank. 
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Appendix  B.  Derivation  of  Equation  50 


Combining  equation  48  from  the  body  of  this  report  and  equations  A-3  and  A-4  from  appendix 
A,  we  get 

—  = — ^1*^2  I  J_/  /  +J_/  /  I 


rs  ^  kl  ^\k^U  ^2khl 


KiX;.,,p':*Bj‘Apyj,,+p)+ 

BJil't  (p”k,l2,+p)pBjAp"‘<ph,Pp)pBjAp~‘ip‘i,  +P)- 


_  D  ^2 

^11  - - 


Uy  +  aj 


I  A  I  n  ^1 

'l2  -°22 

aj  +  aj  <2;  +  aj 


A  straightforward  ealeulation  then  gives  us  the  following; 


j_LAl,(‘Ojv 

271  1-u  " 


2(l-i^,)«2+(l-2vJaia2+2v^a'-(l-2v/Jafa2+  ^ 

a,a2(a,+a2) 

2  (l  -  K J  gf  +  (l  -  2v J  a2^f  +  "(1-2^)  ^2^1  +  /  ^  n 

g,g2(gi+g2)  ^ 


Combining  equation  39  from  the  body  of  this  report  with  equation  B-2,  we  arrive  at  the 
following  equation: 


V  l-2v„ 

—  =  1  +  <■  +  ^ iV,  =  1  +  <• 

K„  A-'’Ap  ' 


2  (1  -  ■'d  <  +  (1  -  2^ )  a.g;  +  2r,a,’  -  (1  -  )  gf g, 

Ijua^a^  i^a^+a^) 

2  (1  -  ^ )  ^1  +  (1  -  2v J  a2^f  +  2v,a2  -  (1  -  2^ )  ^2^1 

Ijua^a^  (gj  +«2) 


P'JLL+P  — 


+  (  ) 
-  p^'kpL+p  —  • 


(B-l) 


(B-2) 


(B-3) 
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Equation  B-3  can  be  rewritten  as  follows; 


+ 


+ 


1 

(N 

\a\+\ 

(i-^J 

\a^al+2v^al  -| 

1  af  aj 

Ijua^ 

(flj  +a2) 

2(1-^) 

\al  +1 

1  a^al  +  2v^al  - 1 

\alay 

Ijua^ 

(<3i  +<32) 

p:%;2,+p 


V^oy 

V^oy 


(B-4) 
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COPIES  ORGANIZATION 


1  MCGILL  UNIVERSITY 

MECHANICAL  ENGINEERING  DEPT 
D  L  FROST 

817  SHERBROOKE  ST  W 
MONTREAL  QUEBEC  H3A  2K6 
CANADA 

1  UNIVERSITY  OF  COIMBRA 
ADAIASSOC 

DEVELOP  INDUSTR  AERODYNAMICS 
I  E  PLAKSIN 

COIMBRA  PORTUGAL  3030-788 

1  UNIVERSITY  OF  LEEDS 
SCHOOL  OF  MECHANICAL 
ENGINEERING 
G  J  SHARPE 
WOODHOUSE  LANE 
LEEDS  WEST  YORKSHIRE  LS2  9JT 
UK 

1  ATOMIC  WEAPONS 
ESTABLISHMENT  UK 
AWE  PEC  ALDERMASTON 
N  K  BOURNE 
READING  RG7  4PR 
UK 

1  ATOMIC  WEAPONS 
ESTABLISHMENT  UK 
AWE  PEC  ALDERMASTON 
J  C  F  MILLET 
READING  RG7  4PR 
UK 

1  UNIVERSITY  OF  CAMBRIDGE 
CAVENDISH  LABORATORY 
C  H  BRAITHWAITE 
CAMBRIDGE  CB3  OHE 
UK 

I  UNIVERSITY  OF  CAMBRIDGE 
CAVENDISH  LABORATORY 
A  L  COLLINS 
CAMBRIDGE  CB3  OHE 
UK 

1  UNIVERSITY  OF  CAMBRIDGE 
CAVENDISH  LABORATORY 
D  J  CHAPMAN 
CAMBRIDGE  CB3  OHE 
UK 


I  UNIVERSITY  OF  CAMBRIDGE 
CAVENDISH  LABORATORY 
D  DRODGE 
CAMBRIDGE  CB3  OHE 
UK 

1  QINETIQ 
M  D  COOK 
BLDG  X59 
FORT  HALSTEAD 
SEVENOAKS  KENT 
TNI 6  7BP 
UK 

I  QINETIQ 
T  A  VINE 
BLDG  X59 
FORT  HALSTEAD 
SEVENOAKS  KENT 
TNI 6  7BP 
UK 

I  SWEDISH  DEFENCE  RESEARCH 
AGENCY 

FOI-SWEDISH  DEFENCE 
RESEARCH  AGENCY 
FOI  GRINDSJON  SE  147  25 
A  HELTE 
TUMBA  14725 
SWEDEN 
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